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DOUBLE ELECTROSTATIC PROBE FOR MEASURING DENSITY, 
TEMPERATURE, AND VELOCITY OF A FLOWING PLASMA 
by Donald L. Chubb 
Lewis Research Center 

SUMMARY 

A method for obtaining plasma electron temperature and density, as well as the 
Mach number and flow velocity from the current- voltage characteristic of a flat-faced 
double probe, is presented. Calculated momentum fluxes of a flowing argon plasma 
obtained with this probe are compared with experimentally determined momentum fluxes. 
Reasonable agreement is obtained. 


INTRODUCTION 

Electrostatic probes have been used to measure plasma properties for many years 
(ref. 1). Double probes and Langmuir probes are used to measure the electron temper- 
ature and plasma number density in stationary plasmas. However, little work has been 
done in high velocity plasmas. In this study a double probe was used to measure the 
plasma Mach number and velocity, as well as the electron temperature and plasma num- 
ber density. 

The following section is an analysis of the double probe. The various plasma prop- 
erties are deduced from the current- voltage characteristic of the double probe. In the 
last section experimental results obtained with the double probe are presented. 

In the analyses it is assumed that the Debye length of the plasma is less than the 
probe characteristic dimension. Also, all collisional processes and magnetic -field 
effects are neglected. Therefore, all collisional mean free paths and the cyclotron radii 
of the electrons and ions must be large compared with the Debye length. 

DETERMINATION OF PLASMA PROPERTIES FROM PROBE DATA 

Figure 1(a) is a schematic diagram of the double probe being considered. One face 



of the probe is alined perpendicular to the plasma flow, and the other face parallel to 
the flow. If a potential is applied by the power supply arrangement shown in figure 1(a), 
the probe current- voltage characteristic will look like figure 1(b). 

Since the whole probe system is "floating, " no net current can flow from the plas- 
ma. For the plasmas to be considered, the ions are "cold" compared with the elec- 
trons. As a result, the maximum current that can flow in the system will be determin- 
ed by the ion saturation current for each face rather than the electron saturation current. 
Since the ion saturation current is much less than the electron saturation current, the 
bulk plasma is much less disturbed by a double probe than a probe that collects electron 
saturation current. Since face 2 is perpendicular to the flow, it will have a larger ion 
saturation current than face 1. As a result, the maximum positive current Ig + is 
greater than the maximum negative current I^ + . In the ideal case Ig + and I^ + are 
constant. In reality the sheath about each face continues to grow, thus giving a larger 
effective current collecting area. As a result the saturation ion currents continue to 
increase slightly as | vJ is increased. 

If there were no flow velocity (U = 0) and both faces of the probe were of equal area, 
the current-voltage characteristic would be symmetric as shown by the dashed curve in 
figure 1(b). The case for U = 0 is discussed by Chen in reference 1 (p. 178-183). 

Now consider how the plasma properties can be deduced from the current-voltage 
characteristic of the probe. Figure 2 is a schematic diagram of the potential distribu- 
tion in the plasma between the two faces of the probe. The potentials Vj and VT, are 
measured with respect to ground and and V2 with respect to the local plasma po- 
tential, Vpj at each probe face. The difference in plasma potential between face 1 and 
face 2 is V c> As a result, 



(Symbols are defined in the appendix. ) 
faces of the probe, 


- V* = V 2 - Vj + V c 

= V pl, 1 " V pZ,2 
Since the same current 


(la) 

(lb) 

must flow through both 


1 ~ i 2+~ i 2- 

(2a) 

1 = H- ~ i l + 

(2b) 


where i 1+ and ig + are the ion currents collected by faces 1 and 2, respectively, and 
i^_ and ig_ are the electron currents collected by faces 1 and 2, respectively. (Eq. (2) 


2 



assumes that the positive direction for current flow is from face 2 to face 1. ) 

Also, since the whole probe system is floating, no net current is collected from the 
plasma, thus 


l l + + x 2 + 


- 1 


1- 


- lr 


= 0 


(3) 


Equation (3) shows that 


x l+ + l 2+ 


= l 


1 - + l 2- 


(4) 


We now assume that, in the sheath surrounding each face, the electron density 
obeys the Boltzmann density distribution (ref. 1, p. 135). Furthermore, we assume 
that the densities outside the sheath on each probe face are the same and that the elec- 
tron temperatures at each face are the same. Thus, 


-qV t AT e 

n l- - n co e 


(5a) 


= n oo e 


-qV 2 AT { 


(5b) 


where n M is the plasma density outside the sheath and T g is the electron tempera- 
ture. The electron current collected by each face is, therefore, 


‘i- = Vi- 


_e 

4 


(6a) 


* 2 - 


= A 2 n 2- V 


(6b) 


where and A 2 are the probe surface areas and v g is the average electron thermal 
speed, which may be expressed as 



Because we are assuming that the Debye length is small compared with the probe radius, 
the collection area for current is just the surface area of the probe face. 
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Using equations (5) and (7) in equation (6) gives 


h- = V- 


kT e -qV 1 AT £ 


2vm, 


(8a) 


i 2- = A 2 n °° 


kT e -qV 2 AT £ 


2nm„ 


(8b) 


In using equation (8) for the electron current collected by the probe, we have ne- 
glected the contribution to i 2 _ from the directed velocity U. However, since 


U « 



for most plasmas this is a good assumption. 

Now substituting equation (8) into equation (4) yields 


H+ + *2+ = 


kT e /. -4V kT e . -qVjAT, 


2jrm, 


A.e 


+ A 0 e 


( 9 ) 


And dividing equation (9) by equation (8b) using equation (la) results in 


i l+ + i 2+ A 1 


QV d qV c , 


2 - 


exp 

A 2 \kT g kT 


+ 1 


( 10 ) 


Substituting equation (2a) for yields 


i l + + i 2 + , A 1 / qV d ^ V c 


- 1 = — exp 


l 2+ - I 


kT kT 

e e; 


( 11 ) 


This is the expression given by Johnson and Matter (ref. 2). 
Taking the logarithm of equation (11) yields 
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(12) 


In 


fl l + + x 2+ 


2 + 



qV c qV d 
— _ + — _ 

kT kT 
e e 


Equation (12) shows that a plot of the left-hand side as a function of should yield a 
straight line. The slope of the line is equal to the inverse of the electron temperature 
in electron volts. In many cases it is assumed that ig + and ij + are constant. How- 
ever, Johnson and Malter (ref. 2) found that ig + and i^ + could be approximated by 
linear functions of V^. Making this approximation, equation (12) can be written as 


where 


AV, + B 

In/ £ 1) = c, + c 9 V. 

^2 V d + b 2-‘ ' 


ii + =a i v d + b i 


f2+ ~ a 2 V d + b 2 


A ~ 3-j + ^2 


B = bj + bg 


A 1 qV c A 1 

Cl = In — - = in — - c 9 V„ 

A 2 kT e A 2 


c 2 


kT, 


(13) 

(14a) 

'j 

(14b) 

(15a) 

(15b) 

(16a) 

(16b) 


The constants a^, bj, ag, and bg can be obtained by fitting a straight line to the ion 
saturation regions of the probe trace. Using the experimental data to be presented later, 
a least squares curve fit (ref. 3) of equation (13) was made to obtain Cj and Cg. In 
this case 


c l = 


X 2 P 1 - X 1 P 2 


(17) 


1^2 - Xj 
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( 18 ) 


where 



(19a) 


(19b) 


(19c) 


(19d) 


The quantity N is the total number of data points. 

A value for the flow Mach number can be obtained from the ratio of front face satu- 
ration ion current to side face saturation ion current. In reference 4, the following re- 
sult for the saturation ion current to a probe in a plasma moving with velocity U is 
presented: 


2 8kT i 

t = V qi '» u + “ 

' Tria- 
ls 


, 1/2 


U > 0 


(20a) 


If U = 0, the Bohm expression for saturation ion current (ref. 1, p. 150) replaces equa- 
tion (20a): 


6 


1 

2 


Apr^oo 



U = 0 


(20b) 


where A pr is the probe surface area and m i is the ion mass. Using equations (20a) 
and (20b) results in 



( 21 ) 


Equation (21) assumes that 



In most cases U 


2 


>> 8kT^/77mp so that 



The definition of Mach number is 




(23) 


where p^ is the plasma density, p^ is the plasma pressure, and y is the ratio of 
specific heats. Since neutrality exists in the plasma, 


Poo = (m e + m t ) « m i n o0 (24) 

Poo = neo (kT i + kT e ) (25a). 

In most cases the ion temperature T i is much less than the electron temperature T g . 

As a result, 
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(25b) 


p ft; n kT 

^oo oo g 

Now substitute equations (24) and (25b) into equation (23). 



m i U 2 


(26) 


Using equation (22) in equation (26) yields 


M 


OO 



(27) 


Equation (27) gives the plasma Mach number in terms of the ion saturation currents. It 
is not necessary to know T g in order to determine the Mach number. A simple meas- 
urement of the saturation ion currents collected by the two faces of the probe is all that 
is necessary. 

With M^ and T g known, the plasma flow velocity U can be calculated from equa- 
tion (26) 


U=M 0C 



(28) 


Also, knowing the saturation ion current I^ + and the electron temperature T , equa- 
tion (20b) can be used to find the plasma density n^ . From (20b), remembering that 
U = 0 for face 1, k 


n 


OO 



(29) 


Equations (16a) and (27) to (29) enable one to calculate the electron temperature, 
plasma Mach number, plasma velocity, and plasma number density from the double 
probe current- voltage characteristic. In the next section these properties will be calcu- 
lated from an experimental probe trace. 
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EXPERIMENTAL RESULTS 


Construction of the flat-faced double probe used in the experiment is illustrated in 
figure 3. Also shown is the electrical circuit. Shields were placed around each face of 
the probe. The shields are biased independently to the same potential as the face they 
enclose. This was done to reduce the amount of current arriving at a probe face from 
any direction other than that perpendicular to the probe face. In other words, to make 
the probe as one-dimensional as possible. The independent biasing of the shields and 
probe faces is accomplished by mounting matched potentiometers on the same shaft. 

The following procedure was used to construct the double probe. A 0. 0508- 
centimeter-diameter tantalum wire was placed inside the 0. 203-centimeter-outside- 
diameter tantalum tube. The space between the wire and tube was then packed with mag- 
nesium oxide. To insure a tight fit between the wire and tube, the tantalum tube was 
then swagged down to an approximately 0. 165 -centimeter outside-diameter. The re- 
sulting tube-wire combination was bent to a 45° angle at one end. The end was then 
machined at a 45° angle to produce one section of the probe. A similar procedure was 
carried out with another wire-tube combination to produce the other section of the probe. 

Each section was then placed in an oval-shaped, two-holed ceramic tube to complete the 

-3 

construction. The area of each face was 2. 87x10 square centimeter. 

/ 

With a magnetoplasmadynamic (MPD) arc (ref. 5) as the source for a flowing argon 
plasma, a typical current-voltage characteristic for the double probe is shown in fig- 
ure 4. The current in the ion saturation regions of the curve increases faster with no 
bias potential on the shields than it does with bias potential on the shields. Biasing the 
shields causes the current to be almost constant in the ion saturation regions. This is 
the expected result for a purely one- dimensional probe. In figure 4 only probe face 1, 
which is parallel to the flow, shows the effect of a biased shield. However, in other 
cases the effect has been noted in both ion saturation regions. The probe curve was 
unchanged outside the ion saturation regions with or without shield biasing. The plasma 
parameters are determined mainly by the region of the probe curve between the ion 
saturation regions. Therefore, we conclude that any potential difference that may de- 
velop between one of the probe faces and its corresponding shield as a result of the in- 
dependent biasing of the shields has negligible effect on the plasma parameters deduced 
from the probe curve. 

Probe current-voltage data were sent directly to an IBM 360 computer. Program- . 
ming equations (16) to (19) and (27) to (29) allowed on-line computation of the electron 
temperature, plasma Mach number, velocity, and density. The saturation currents, Ig + 
and I 1+ , were taken at the points where the linear saturation portions of the probe trace 
(fig. 4) begin to deviate from the exponentially varying middle portion of the trace. Fig- 
ure 5 compares the experimentally determined left-hand side of equation (13) with the 
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least squares approximation of this function given by the right-hand side of equation (13). 
Slightly different results occur depending on the number of data points used in the least 
squares approximation. Including more points in the ion saturation portion of the 
current-voltage characteristic produced higher electron temperatures and a poorer fit 
in the region between the ion saturation portions. Including points for -6 s V d s 12 
yielded kT g /q = 2. 01 electron volts, and points for -2.5 ^ 7.5 yielded kT e /q - 

1. 80 electron volts. Table I summarizes the results obtained from the data given in 
figure 4. 

To check the double probe results, a comparison was made between the momentum 
o 

flux, m-mU , calculated from double probe results, and the momentum flux, deter- 
mined by measuring the thrust of the MPD arc. Because the Mach number of the flow 
is high, the pressure term is small when compared with the momentum flux term in 
the equation of motion. Therefore, the thrust produced by the MPD arc is just the inte- 
grated momentum flux in the beam: 

T = J m.nv 2 dA, + f m n v 2 dA, (30) 

*'*1 b o o o b 

A b A b 


The first integral is the thrust contribution from the ions, and the second integral the 
contribution from the neutrals. Since the electron mass is small, the electron momen- 
tum has been neglected in equation (30). Terms appearing in equation (30) are the ion 
mass mj, the electron density n, the ion velocity v, the neutral mass m Q , the neutral 
density n Q , the neutral velocity v Q , and the plasma beam cross sectional area A^. In 
using the electron density in equation (30) we are assuming charge neutrality. 

If we assume that the density and velocity are independent of azimuthal angle and 
that « m Q , then equation (30) becomes 

2 2 

T = 277m. I nv r dr + 277m. I n.v^r dr (31) 

\Jq i •'o o o 


where R is the radius of the plasma beam . In order to carry out the integrations in 
equation (31), the radial dependence of the densities and velocities must be known. The 
double probe can be used to determine the ion velocity and density distributions. How- 
ever, there is no convenient way to obtain the neutral density and velocity distributions. 
As a result, we approximated the neutral thrust contribution as follows: 


T o 



n v 2 r dr 
o o 


m U 
o o 


(32) 
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where m Q is the neutral mass flow rate and U 0 is some average neutral velocity. 
Using equation (32) in (31) yields the following: 

fR n 

T = 2nm. / mTr dr + m U (33) 

1 Jq O O 

The neutral flow rate can be obtained from the measured total mass flow rate m t and 
the calculated ion flow rate m^: 

r R 

m Q = m t - m i = m^ - 2vm^ Jq nvr dr (34) 

Substituting equation (34) into (33) yields 

f R 2 

T = 2jrm- / nv r dr + 

X Jq O 

Double probe data for the density and velocity distributions at an axial station 
61 centimeters from the exit of the MPD arc is shown in figures 6 and 7. We approxi- 
mated these data with the following functions: 


- 2?rm 


f 

i 


nvr dr 


(35) 


n=n TO e- 

(36) 

v = u(i - 

(37) 

\ R/ 



where n M and U are the density and velocity in the center of the beam. Fitting the 
data in figures 6 and 7 with the straight lines shown, we find that R « 40. 7 centimeters 
and 1/a « 6.27 centimeters. 

Substituting equations (36) and (37) into (35) yields the following result: 



For argon with 1/a = 6. 27 centimeters and R = 40. 7 centimeters, equation (38) be- 
comes 
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(39) 


T - 8.53X10 -28 n M U 2 + U Q (m t - 1.64X10" 27 n M U) (SI units) 

The ion velocity U and density n M are obtained from the double probe results. How- 
ever, we have no measurement of the neutral average velocity U Q . Certainly U Q ^ U, 
since the neutrals are accelerated by collisions with the ions. Sovie and Connolly 
(ref. 6) measured ion and neutral velocities in the exhaust of an MPD arc running on 
ammonia. Based on their results, we approximated the average neutral velocity as 
U Q « 2/3 U. Using U Q = 2/3 U as well as the velocities and densities given in table I 
in equation (39) yields the results for thrust shown in tablell. The experimentally 
determined thrust was 0. 28 newton. 

A comparison of the calculated and measured thrusts shows that there is reasonable 
agreement between them. The calculated thrust is about 30 percent less than the meas- 
ured thrust. If we assume that U Q « U then the measured and calculated thrusts are 
almost equal. Therefore, within the accuracy of the neutral thrust approximation, we 
conclude that the plasma properties determined by the double probe are a good approxi- 
mation to actual plasma properties. 


CONCLUSIONS 

A method is presented for calculating the plasma electron temperature, number 
density, flow Mach number, and flow velocity based on the current-voltage character- 
istic of a flat-faced double probe. The analysis assumes a Maxwellian distribution func- 
tion for the electrons. Experimental results for the plasma properties using a double 
probe were obtained in the exhaust of an magnetoplasmadynamic (MPD) arc. The thrust 
calculated using the double probe determined plasma properties was compared with the 
measured thrust. Reasonable agreement was found, indicating that the double-probe 
determined plasma properties are a good approximation to the actual plasma properties. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 17, 1973, 

503-10. 
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APPENDIX - SYMBOLS 

A constant defined by eq. (15a) 

cross sectional area of MPD arc beam 
Apr probe surface area 
A sheath area 

Aj surface area of face 1 of double probe 

Ag surface area of face 1 of double probe 

B constant defined by eq. (15b) 

c ^ constant defined by eqs. (16a) and (17) 

c 2 constant defined by ieqs. (16b) and (18) 

I probe current 

Ij + saturation ion current to probe face 1 

I 2+ saturation ion current to probe face 2 

i 1+ ion current to probe face 1 

i 2+ ion current to probe face 2 

ij_ electron current to probe face 1 
i 2 _ electron current to probe face 2 
i_ electron current to Langmuir probe 
i + ion current to Langmuir probe 
k Boltzmann constant (1.38X10 -22 J/K) 

M m Mach number outside probe sheath 
m mass 

m mass flow rate 

n^ electron number density in sheath about probe face 1 

n 2 electron number density in sheath about probe face 2 

n w electron number density outside probe sheath 
Pj defined by eq. (19c) 

P 2 defined by eq. (19d) 

p pressure 



q magnitude of electron charge (1. 602 C) 
r radius in cylindrical coordinate system 

T total thrust 

T cal total calculated thrust 

T g electron temperature 

T q neutral thrust 

U velocity of ion flow 

U Q velocity of neutral flow 

V electric potential 

V c difference in plasma potential between probe face 1 and face 2 
V d difference in potential between probe face 2 and face 1 of double probe 
V p £ plasma potential 

V I potential of probe face 1 with respect to plasma potential 

V 2 potential of probe face 2 with respect to plasma potential 

V'^ potential of probe face 1 with respect to ground potential 

Vr> potential of probe face 2 with respect to ground potential 

v velocity 

v average electron thermal speed 

Xj defined by eq. (19a) 

X 2 defined by eq. (19b) 

y ratio of specific heats 

e Q permittivity of free space (8. 855x10" F/m) 
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TABLE I. - COMPARISON OF PLASMA PROPERTIES DETERMINED FROM 


DOUBLE- PROBE CURRENT- VOLTAGE CHARACTERISTICS 


[Data obtained from fig. 4.] 


Range of data 

Saturation ion current, mA 

Electron 

Mach 

Ion flow 

Number 

points con- 

J l + 

J 2 + 

temperature, 

number, 

velocity, 

density, 

sidered, 

V 

T 

e’ 

eV 

M„ 

U, 

m/sec 

n ~-3 

m 

-6 £ V d s 12 

0.040 

0. 490 

2. 01 

4.72 

1.34xl0 4 

7. 99X10 17 

-2.5 s V d £ 7.5 

.040 

.490 

1. 80 

4.72 

1.27X10 4 

8. 42xl0 17 


TABLE n. - COMPARISON OF THRUST CALCULATED FROM 
DOUBLE PROBE RESULTS AND MEASURED THRUST 


[ion mass flow rate, 1.76x10 3 kg/sec; neutral mass flow rate, 
_ 0 

9x10 kg/sec; total mass flow rate, 2. 66 kg/sec; measured 
thrust, 0.28 n!] 


Number 

Ion 

Neutral 

Ion 

Neutral 

Total thrust, 

density, 

velocity, 

velocity, 

thrust, 

thrust, 

T cal = T i + T o> 

n , 

u, 

U = — U 

T. , 

T , 

N 

°°o 



l* 

o’ 


m 

m/sec 

o 

N 

N 


7.99X10 17 

1.34X10 4 

8.94X10 3 

mm 

B 

0.203 

8.42 

1.27 

8.47 

E9 

m 

. 191 




















Experimental points 

Range Electron 
of V d , temper- 
V ature, 




Viu 'i 'luajjno aqoJc) 


-8 

o. 

.-S' 


■s 


T3 

IS 
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Difference in potential between probe 
face 1 and face 2 . V d , V 

Figure 5. - Comparison of experimental data and 
calculated electron temperatures. (Data from 




0 4 8 12 16 20 

Radial coordinate, r, cm 


Figure 6. - Radial density profile in magneto- 
plasmadynamic arc. (nln^ = e' ar . ) 



Radial coordinate, r, cm 

Figure 7. - Radial velocity profile in magneto- 
plasmadynamic arc. ^ ■ 1 - 
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